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Momentum-space azimuthal harmonic event planes (EP) are constructed from final-state mid¬ 
rapidity particles binned in transverse momentum (pr) in ^fSNN = 200 GeV Au-|-Au collisions in 
a multi-phase transport (AMPT) model. The EP correlations between pt bins, corrected by EP 
resolutions, are smaller than unity. This indicates that the EP’s decorrelate over pT in AMPT, 
qualitatively consistent with data and hydrodynamic calculations. It is further found that the EP 
correlations approximately factorize into single pT-bin EP correlations to a common plane. This 
common plane appears to be the momentum-space EP integrated over all pr, not the configuration- 
space participant plane (PP). 

PACS numbers: 25.75.Ld, 25.75.Dw, 24.10.Jv, 24.10.Lx 


I. INTRODUCTION 

In high energy heavy ion collisions, the initial high en¬ 
ergy density and pressure buildup drive the collision sys¬ 
tem to rapid expansion. In non-central collisions, the 
pressure gradient is anisotropic due to the anisotropic 
transverse overlap geometry, resulting in an observable 
momentum-space azimuthal anisotropy [I[. The particle 
azimuthal distribution can be expressed in Fourier series: 

oc 1 + ^2u„cos[n((^ - T'n)], (1) 

where 4> represents the particle azimuthal angle and is 
the harmonic plane angle. The Fourier coefficients 0 

= (cos[n((^-T'n)]), (2) 

characterize the magnitudes of the harmonic 
anisotropies, where the bracket (...) denotes aver¬ 
ages over particles and events. In the limit of smooth 
nuclei, odd harmonics vanish; A ’2 coincides with the 
reaction plane (the plane defined by the beam direction 
and the impact parameter vector), and elliptic flow (^ 2 ) 
is the leading component because of the predominant 
elliptic overlap shape. Large V 2 has been measured 
experimentally Q. In fact, the measured V 2 is so large 
that hydrodynamic descriptions are applicable, and a 
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very small shear viscosity to entropy density ratio (rj/s) 
is required [3|. This constitutes an important evidence 
for the formation of a new form of matter, the strongly 
interacting quark-gluon plasma (sQGP) in relativistic 
heavy ion collisions Q. Due to fluctuations, the energy 
density in the overlap region is lumpy, so all harmonics 
can exist Q and the harmonic planes of different orders 
are unnecessarily the same. Higher order harmonics 
are more sensitive to 77 /s, and hence play an important 
role in extracting precise information about the sQGP 
medium Q. 

As the initial configuration space information of the 
overlap region is not experimentally accessible, the har¬ 
monic symmetry plane in configuration space, 

called the participant plane (PP), cannot be measured. 
In experiment, anisotropic flow is measured via final- 
state particle correlations [^. For example, one con¬ 
structs an event plane (EP), by the final-state par¬ 

ticle momenta, as a proxy for 'I'^^ [7|. One then corre¬ 
lates a test particle with the to measure anisotropic 
flow via Eq.(2) where the deviation of from T'n is 
corrected by a resolution factor Q. Alternatively, the 
anisotropic flow is measured using azimuthal correlations 
between the observed particles, where the Fourier coeffi¬ 
cients of the correlation function are given by 

Vu{pt,Pt) = (cos[7r(^“ - /)]). (3) 

Often factorization is assumed where the two-particle 
correlations arise only from the single particle correla¬ 
tions to a common harmonic plane Tn (i.e. nonflow in¬ 
trinsic particle correlations are neglected). Under such 
conditions, since the common A’n cancel in ((/“ - Ikn) - 
- '!'„), one has 
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^niPTTPx) ~ 'Vn{PT)'^ri{PT)- ( 4 ) 

The single-particle anisotropic flow can be simply ob¬ 
tained as the square root of the two-particle cumulant 
measurement in Eq.( 3 ) by choosing the two particles from 
the same phase space, or when the Vn of one particle is 
known (referred to as the reference particle). 

Studies have shown, however, that the harmonic 
planes, whether in configuration or momentum space, 
are not the same at different pseudorapidities. Previ¬ 
ous works @ have shown that EP’s are decorrelated 
over pseudorapidity (77) and the decorrelation increases 
with increasing 77 gap between the two particles. 

Not only does the harmonic plane depend on 77, but 
also on transverse momentum (pr)- Gardim et al. 0 
have pointed out with ideal hydrodynamic calculations 
that the fluctuating flow angles depend on px- Vis¬ 
cous hydrodynamic calculations m indicate that the 
harmonic planes decorrelate over px even in the same 
77 region. This is indeed observed by CMS in Pb-fPb 
and p-fPb collisions at the LHC [l^- The experimental 
confirmation of the px decorrelation is taken as a strong 
evidence supporting hydrodynamic descriptions of rela¬ 
tivistic heavy ion collisions. 

Experimentally, in order to reduce nonflow correlations 
(such as small angle jet correlations and resonance de¬ 
cays), one often apply 77 gap between the two particles 
used in correlation measurements. In measuring high px 
particle anisotropy, the reference particle is often taken 
at low Px- The harmonic plane decorrelation over 77 and 
Px casts a problem in those measurements; it may require 
the EP to be constructed in the same phase space as the 
particle of interest in order to measure flow anisotropy Q . 
As an operational definition, one may define the flow of 
particles in a particular phase space as that with respect 
to the harmonic plane of the same phase space region. 
Namely, 

Vniv, Pt) = {cos[n{(j>{T],px) - 'Sn{V,PT))])- ( 5 ) 

We note, however, that such Vn is, in realty, contami¬ 
nated by non-negligible nonflow effects. 

The 77 and px decorrelation of the harmonic planes 
explicitly breaks factorization of Eq.( 4 ). This is easy to 
see. The two-particle cumulant measurement with the 
two particles coming from different phase spaces is 

VniPx^Pr) = {cos[n{(l)nipx) - ^n{Px))]) 

= VniPx)VniPx) n{Px) “ ^^(Pt))]- 

( 6 ) 

The two harmonic planes from different phase spaces are 
now different and can no longer cancel in two-particle 
azimuthal correlations. Note that in Eq.(6) we have ne¬ 
glected possible nonflow contributions (which is usually 


small when a large phase-space gap is imposed between 
the two particles), so that we can write the average of 
the product into the product of averages. 

The EP decorrelations were experimentally measured 
to be appreciable [l^ . However, the experimentally mea¬ 
sured two-particle cumulants obey factorization exceed¬ 
ingly well [ 1 ^. It implies that cos[?r(d'„(p!^) — d'„(p^))] 
correlations between two phase spaces, a and b, must fac¬ 
torize to a good degree [ 1 ^. This is by no means obvious. 

In this paper we try to address two questions. One, are 
EP’s in parton transport model similarly decorrelated as 
in data and hydrodynamic calculations? Second, how 
well are the EP correlations factorized? The answers to 
these questions can hopefully reveal additional insights 
to the physics of EP decorrelations. 

II. ANALYSIS METHOD 

The AMPT (A Multi-Phase Transport) model with 
string melting describes many experimental data rea¬ 
sonably well, particularly the anisotropic flow measure¬ 
ments [1^. We thus use the AMPT model with string 
melting for our study. There are four main compo¬ 
nents in AMPT: the initial conditions, parton interac¬ 
tions, hadronization, and hadron interactions. The ini¬ 
tial conditions are obtained from the HIJING model [l^ . 
which includes the spatial and momentum information of 
minijet partons from hard processes and strings from soft 
processes. The time evolution of partons is then treated 
according to the ZPC parton cascade model 0 . After 
parton interactions cease, a combined coalescence and 
string fragmentation model is used for the hadronization 
of partons. The scattering among the resulting hadrons 
is described by a relativistic transport (ART) model [l^ 
which includes baryon-baryon, baryon-meson and meson- 
meson elastic, and inelastic scatterings. 

In AMPT model, two harmonic planes can be con¬ 
structed. One is calculated from the initial configuration 
space information of partons by 

^pp ^ atan2((r"sin(n((ii.)), (r" cos(n(()i.))) -f tt 

where r and are the polar coordinate position of each 
parton. Because of event-by-event geometry fluctuations, 
47^^ is not necessarily as same as the reaction plane. 
Experimentally, the coordinate space information is not 
accessible. Instead, the EP can be constructed from hnal- 
state particle momenta by 

^EP ^ atan 2 ((sin(n(())),(cos(n(()))) 

where (j) is the azimuthal angle of the particle momen¬ 
tum. Due to the finite multiplicity of constituents, the 
constructed harmonic plane is smeared from the true one, 
-the geometry harmonic plane of the participant 
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Pt (GeV/c) 

n ^ 2 

n — 3 

Pt (GeV/c) 

n ^2 

n — 3 

0 - 0.2 

0.11 ± 0.03 

0.06 ± 0.04 

0.2 - 2.0 

0.889 ± 0.001 

0.568 ± 0.004 

0.2 - 0.4 

0.557 ±0.005 

0.17 ± 0.02 

0 - 0.2 II 0.4 - 2.0 

0.853 ± 0.002 

0.537 ± 0.004 

0.4 - 0.6 

0.719 ± 0.003 

0.339 ± 0.008 

0 - 0.4 II 0.6 - 2.0 

0.810 ± 0.002 

0.435 ± 0.006 

0.6 - 0.8 

0.684 ±0.003 

0.384 ± 0.007 

0 - 0.6 II 0.8 - 2.0 

0.829 ± 0.002 

0.435 ± 0.006 

0.8 - 1.0 

0.608 ± 0.004 

0.364 ± 0.007 

0 - 0.8 II 1.0 - 2.0 

0.847 ±0.002 

0.462 ± 0.005 

1.0 - 1.2 

0.525 ± 0.006 

0.313 ± 0.009 

0 - 1.0 II 1.2 - 2.0 

0.858 ± 0.001 

0.480 ± 0.005 

1.2 - 1.4 

0.423 ± 0.008 

0.26 ± 0.01 

0 - 1.2 II 1.4 - 2.0 

0.864 ± 0.001 

0.496 ± 0.005 

1.4 - 1.6 

0.36 ± 0.01 

0.23 ± 0.01 

0 -1.4 II 1.6 - 2.0 

0.866 ± 0.001 

0.506 ± 0.005 

1.6 - 1.8 

0.32 ± 0.01 

0.21 ± 0.01 

0 - 1.6 II 1.8 - 2.0 

0.867 ±0.001 

0.507 ± 0.005 

1.8 - 2.0 

0.33 ± 0.01 

0.23 ± 0.01 

0 - 1.8 

0.867 ±0.001 

0.507 ± 0.005 


TABLE I: (Color online) The resolutions, of the second- and third-order harmonic EP’s constructed by final-state particle 
momentum azimuthal angles in different transverse momentum {pr) bins. Particles are taken from pseudo-rapidity range of 
I 77 I < 1.0 in 20-40% centrality (b = 6.79-9.61 fm) Au-|-Au collisions at ^snn = 200 GeV simulated by the AMPT model (string 
melting, 3 mb parton cross section). 


\ Pt 

Pt 

0 - 0.2 

0.2 - 0.4 

0.4 - 0.6 

0.6 - 0.8 

0.8 - 1.0 

1.0 - 1.2 

1.2 - 1.4 

1.4 - 1.6 

1.6 - 1.8 

1.8 - 2.0 

1.8 - 2.0 

-0.32 ± 0.38 

0.79 ± 0.14 

0.70 ± 0.07 

0.82 ± 0.06 

0.86 ± 0.07 

0.92 ± 0.08 

0.90 ± 0.09 

0.70 ± 0.12 

0.65 ± 0.12 


1.6 - 1.8 

-0.43 ± 0.43 

0.85 ± 0.16 

0.87 ± 0.08 

0.91 ± 0.07 

0.77 ± 0.08 

0.92 ± 0.09 

0.81 ± 0.11 

0.98 ± 0.13 


0.86 ± 0.05 

1.4 - 1.6 

-0.79 ± 0.43 

0.51 ± 0.16 

1.00 ± 0.08 

1.07 ±0.07 

0.98 ± 0.08 

1.10 ± 0.09 

1.06 ± 0.11 


0.94 ±0.05 

0.90 ± 0.05 

1.2 - 1.4 

0.16 ± 0.34 

0.66 ± 0.13 

1.00 ± 0.06 

1.00 ± 0.06 

1.06 ± 0.06 

0.87 ±0.07 


0.97 ± 0.04 

0.93 ± 0.04 

0.92 ± 0.04 

1.0 - 1.2 

-0.32 ± 0.28 

0.88 ± 0.11 

0.85 ± 0.05 

0.88 ± 0.05 

1.03 ± 0.05 


1.02 ± 0.02 

0.95 ± 0.03 

0.88 ± 0.03 

0.85 ± 0.03 

0.8 - 1.0 

-0.32 ± 0.24 

0.68 ± 0.09 

1.00 ± 0.05 

0.88 ± 0.04 


0.98 ± 0.02 

0.95 ± 0.02 

0.97 ± 0.03 

0.93 ± 0.03 

0.86 ± 0.03 

0.6 - 0.8 

-1.00 ± 0.23 

0.82 ± 0.09 

0.88 ± 0.04 


0.95 ± 0.01 

0.97 ±0.01 

0.96 ± 0.02 

0.92 ± 0.02 

0.93 ± 0.03 

0.89 ± 0.02 

0.4 - 0.6 

-0.45 ± 0.26 

0.76 ± 0.10 


0.97 ±0.01 

0.96 ± 0.01 

0.95 ± 0.01 

0.94 ± 0.02 

0.90 ± 0.02 

0.91 ± 0.02 

0.86 ± 0.02 

0.2 - 0.4 

0.60 ± 0.51 


0.95 ± 0.01 

0.94 ±0.01 

0.95 ± 0.02 

0.97 ±0.02 

0.93 ± 0.02 

0.92 ± 0.03 

0.96 ± 0.03 

0.87 ±0.03 

0 - 0.2 


0.81 ± 0.09 

0.97 ± 0.07 

0.87 ±0.08 

0.92 ± 0.09 

1.06 ± 0.10 

1.00 ± 0.12 

0.94 ±0.14 

0.99 ± 0.16 

0.86 ± 0.16 


TABLE IL (Color online) The resolution-corrected EP correlation strengths C2{pt,Pt) (right down corner) and C3(j)t,Pt) 
(left up corner) of Eq.(9) for different combinations of p^ and p^ {p^ 7 ^ Pt)- Particles used in EP construction are restricted 
within pseudo-rapidity |r;| < 1.0. The data are simulated by AMPT (string melting, 3 mb cross section) in 20-40% centrality 
(b = 6.79-9.61 fm) in Au+Au collisions at ^/snn = 200 GeV. 


partons in configuration space in the limit of infinite par- 
ton multiplicity- by a resolution factor. The resolution 
factor Kn = (cos[n(l['n — “'^)]) is calculated with an 

iterative procedure by the sub-event method, dividing 
the constituents randomly into two sub-events Q- Be¬ 
cause of the large initial parton multiplicity, the calcu¬ 
lated initial configuration-space PP resolution, is 

nearly unity However, the final-state hadron multi¬ 
plicity is smaller so the final-state momentum-space EP 


resolution, is smaller than unity. The EP resolu¬ 

tions in different pp bins are tabulated in Table I for both 
the second- and third-order harmonic EP’s. In addition, 
the resolutions of the EP’s constructed from all particles 
in 0 < Pt < 2 GeV/c but excluding those in a particular 
Pt bin are also listed. 
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FIG. 1: (Color online) Upper panel shows the EP correlation strength Cn{pT,P^), corrected by the corresponding EP resolu¬ 
tions, as a function of py for a fixed p^ bin. Particles used in EP construction are restricted within pseudo-rapidity |t;| < 1.0. 
The data are simulated by AMPT (string melting, 3 mb cross section) in 20-40% centrality (b = 6.79-9.61 fm) in Au-|-Au 
collisions at ^snn = 200 GeV. Red lines are the fit to data by the function of CniPr) x CniPr) where Cnipr) is a set of ten fit 
parameters. Lower panel shows the ratio of the data points to the fit result. 


III. RESULTS AND DISCUSSIONS 

The EP angles 4 '®^(pt) are constructed from final- 
state charged particle momenta in ten equal-width pT 
bins in 0 < pt < 2 GeV/c. Particles used in the EP 
construction are restricted within \r]\ < 1.0. Both the 
elliptic (n = 2) and triangular (n = 3) EP’s are con¬ 
structed. The EP correlation between and p^ bins, 
corrected by the EP resolutions in Table I, is given by: 


Cn{PTTP\') 


(cos[n(4^rbT)-^rb^))]) 


( 9 ) 


The EP correlation strengths Cn are tabulated in Ta¬ 
ble II for both n = 2 and n = 3. For each CniPTiPr) 
value in the table,we have already taken the average of 
CniPT^Pr) CniPT^Pr)- from the table,the 

measured CniPTiPr) Pt 7^ Pt 1®®® than unity. 
This indicates that the EP’s in AMPT are indeed decor- 
related over Pt- The decorrelations are similar in mag¬ 
nitude to those observed in data at the LHC [l^ and 
from hydrodynamic calculations [HI- Figure 1(a) shows 
the EP correlation strengths Cn as a function of Pt for 
a fixed p\, = 0.4 < px < 0.6 GeV/c, as an example. 


To test factorizability, the EP correlations of all combi¬ 
nations of p^ and py are fitted by a set of ten parameter 
Cnipr) corresponding to the ten px bins: 


CniPTyPr ) = Cnipx ) X CniPx ) ■ (10) 

Factorization of the EP correlations can be tested by the 
fitting quality. The fit result, c„(py) x Cn(py), is shown 
as the curves in Fig.l along with the EP correlations data 
points from AMPT, CniPxiPx)- l^^l® ®''^® good for both 
n=2 and n=3. To illustrate the fit quality. Fig. 1(b) 
shows the ratio of the AMPT data to the fit. The ratio 
is consistent with unity. Note Fig.l is only a subset of the 
EP correlation data for a fixed p^ bin as an example, but 
all 10 X 9/2 correlation data for a given harmonic (n = 2 
or 3) are fitted simultaneously with ten fit parameters in 
a single fit. The fit y^/ndf are 30.5/35 and 45.5/35 for 
C 2 and C 3 , respectively. Our fit results indicate that the 
EP correlations are well factorized. The EP correlation 
between two px bins is determined by only the properties 
of the two single-pT bins. This may indicate that the EP 
decorrelation is of a random nature. 

Since CniPxiPx) factorize well into products of single 
PT-bin EP quantities, the fit results suggest that the EP 
correlations may be caused by single py-bin EP correla- 
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FIG. 2: (Color online) Upper panels show the comparisons of the fit parameter c„{pt) to Cn(j)T, the correlations between 

'^n^iPr) and the EP constructed from all pr (but excluding the relevant pr bin) and Cn{pT,'i^n^), the correlation between 
'^n^iPr) and the PP, for (a) second- and (b) third-order harmonic. Lower panels show the corresponding ratios. Data are for 
20-40% centrality (b = 6.79-9.61 fm) in Au-|-Au collisions at ^/snn = 200 GeV from the AMPT model (string melting, 3 mb 
parton cross section). 

tions to a common plane, Namely, 


J 


(cos[n(vl/r(p^)-^r(P^))]) 


(cos[n(^g^(pg.) - ^„)]) (cos[n(^g^(p^) - ^„)]) 

^riPr) 


( 11 ) 


and 


c „( pt ) 


(cos[n(^g^(pr) - ^n)]) 
^S^{pt) 


( 12 ) 


The good fitting quality demands an answer to the 
following question: what plane is 'I'n? To address this 


question, we show in upper panel of Fig. 2 the compar¬ 
isons between the fit parameter c„(pt) and that given 
by Eq.(12) where is in one case calculated by the 
configuration-space and in the other case by the 

momentum-space from particles over the entire pp 
range: 



























and 
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Cn{pT,'i>n^) 


(cos[n(vI/r(PT)-vl^D]) 


(13) 


/ ^ (cos[n(^g^(pT) - <PT <2 GeV/c but exclude pt))]) 

Cn Pt, n < pt < 2 GeV/c but exclude pt) 


respectively. The are constructed from the ini¬ 

tial parton configuration within the pseudo-rapidity of 
\Pparton\ < 1-0. The are constructed from the 

final-state particle momenta within \phadron\ < 1-0 and 
0 < pt < 2 GeV/c but excluding the pr bin of inter¬ 
est. To remove auto-correlations, the construction of 
'1'^^ excludes particles of the particular pr bin used for 
the construction of 'I'„(pt). The resolutions and 

5ft®^(0 < pt <2 GeV/c but exclude pr) are to take into 
account the inaccuracy of the constructed PP and EP 
(although Ift®® is found to be approximately unity). 

The comparisons of c„(pt) to Cn{pT,'^n^) and 
Cn{PT, 'k®®) are shown in Fig.2 for both the second- and 
third-order harmonic EP’s. We find that c„(pT,'k®®) 
is more consistent with the fit parameter c„(pr) than 
Cn(PT,'I'®®) is. The more quantitative comparisons are 
shown in Fig. 2(c) and (d) where the ratios are depicted. 
The ratio of cz{pt, '^ 3 ^)/cz{pt) is consistent with unity 
while C 3 (pt, ^'®®)/c 3 (pt) is not. For n = 2, both ra¬ 
tios are smaller than unity, however, C 2 {pt, 'I'f'®)/c 2 (pT) 
is closer to unity than C 2 {pt,'^ 2 ^)/c 2 {pt)- The devia¬ 
tion of C 2 {pt,'^ 2 ^)/c 2 {pt) from unity could be due to 
nonflow effect where some particles are intrinsically cor¬ 
related in the final-state momentum space. Our results 
show that the global harmonic plane is not the 'I'®®; it 
may be the 'k®®. In other words, 'I'®®(pt) fluctuates 
randomly around 'I'®®, not 4'®®. 


IV. CONCLUSIONS 


In summary, we have studied momentum-space event 
plane (EP) correlations in b = 6.79-9.61 fm Au-|-Au colli¬ 
sions at = 200 GeV using the AMPT model (string 

melting, 3 mb parton cross section). The EP’s are con¬ 
structed from final-state mid-rapidity particles binned in 
Pt- Both the second- and third-order harmonic EP’s 
are studied. The EP correlations between p^ bins are 
corrected by the corresponding EP resolutions. The EP 
correlations are found to be smaller than unity, indicat¬ 
ing that the EP’s decorrelate over px in AMPT, quali¬ 
tatively consistent with data and hydrodynamic calcula¬ 
tions. It is further found that the EP correlations fac¬ 
torize into single p^-bin EP correlations to a common 
plane. This common plane appears to be approximately 
the momentum-space EP integrated over all pt, but not 
the configuration-space participant plane (PP). 
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